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Abstract: We present an optimization of the reduction algorithm of one- loop amplitudes in terms 
of master integrals. It is based on the exploitation of the polynomial structure of the integrand 
when evaluated at values of the loop-momentum fulfilling multiple cut-conditions, as emerged in the 
OPP-method. The reconstruction of the polynomials, needed for the complete reduction, is rended 
very versatile by using a projection-technique based on the Discrete Fourier Transform. The novel 
implementation is applied in the context of the NLO QCD corrections to ud W + W~W + . 
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1. Introduction 

In the last few years, we have been witnessing a boost in developing new ideas aiming to the efficient 
computation of one-loop amplitudes H, as extensively reported in [0. Besides standard techniques, 
where the tensor reduction is explicitly performed, new numerical and analytical developments, 
originally inspired by unitarity arguments ||, f§] have emerged. The common features of the so called 
unitarity-based methods -|P6[, is the change of perspective they propose: instead of focusing on 
the actual evaluation of complete integrals, they pursue the determination of the coefficients of 
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the scalar one-loop functions contributing to the unknown integrals. This possibility relies on the 
fact that the basis of scalar function to express any one-loop integrals is known in terms of Boxes, 
Triangles, Bubbles and (in massive theories) Tadpoles [PTfl . Schematically, one can write a Master 
Equation for any one-loop amplitude M. such as: 

M. = di Boxj + Cj Triangle^ + bi Bubble j + a« Tadpole^ , (1.1) 

i i i i 

where di, Ci, bi and a; are the coefficients to be determined. 



Very recently in 28, E{J, by exploiting the properties of spinor- integration of double-cuts of di- 



mensionally regulated integrals |10|, [L2|, [T|, [14], [15], [15|, |17|, general analytic formulas for the 



coefficients di, Ci, bi were presented. These formulas can be evaluated at the occurrence, without 
performing any integration, by specializing the value of input variables that are specific to the initial 
cut-integrand, which is assembled from tree-level amplitudes. 

Alternatively to any phase-space integration, in [JTSL |19[ it was proposed a very efficient method 
for the reconstruction of the so-called 4-dimensional cut-constructible term of any scattering am- 
plitude, corresponding to the poly- logarithmic structure arising when Eq. ( |1.1|) is expanded around 
4-dimensions. This method, by-now known as OPP-reduction, allows the numerical reconstruction 
of (the 4-dimensional limit of) the coefficients, di, Ci, bi and <2j, by solving a system of algebraic 
equations that are obtained by: i) the numerical evaluation of the integrand at explicit values of the 
loop- variable, on the one side; ii) and the knowledge of the most general polynomial structure of the 



integrand itself ||30|| , on the other one. We remark that the values of the loop momentum used for 
the numerical evaluation of the integrand are chosen among the set of solutions of the multiple-cut 
conditions, namely the solutions of the system of equations obtained by imposing the vanishing of 
the denominators on each 4-dimensional cut. 

For the complete evaluation of scattering amplitudes, one has to consider that the 4-dimensional 
expansion of Eq. ( |1 . 1| ) generates not only a poly- logarithmic term, but as well a rational term 



which cannot be detected by (massless) cuts in 4-dimensions. In [21], it has been recently shown 
that there are two sources of the rational terms: the first contribution, that is quite simple to 
calculate, originates from the generic (n — 4)-dimensional structure of the numerator of any one-loop 
amplitude and it can be derived by using appropriate Feynman rules within a tree-like computation. 
The second contribution originates instead from the reduction of the 4— dimensional part of the 
numerator in terms of the n— dimensional denominators appearing in the scalar integrals. This 
part, that is more subtle to extract, can be computed within the OPP-method in a completely 



automatized way, following one of the approaches discussed in [^IJ and numerically implemented in 
the public code CutTools 



Alternatively, the reconstruction of the rational term can be achieved as well by using techniques 
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like direct computation [31, [32|], by the bootstrapping method by cuts in n-dimensions [23|, p9| . 



or by explicitly computing the amplitude at different integer value of the space-time dimensions [^6 
In particular, in |25|, |26| it has been proposed an extension of the OPP-reduction, implementing an 
integrand decomposition valid in n-dimension, rather than in 4-dimension, which exposed a richer, 
yet polynomial, structure of the cut-integrand. 

The efficiency of the OPP-reduction has been shown in non-trivial applications, like the 6-photon 
amplitudes with massless and massive fermion-loop ||19|| , the virtual QCD correction to qq — > ZZZ 
0], and the complete cross section for the production of three vector bosons at LHC pi 



Within the OPP-reduction, the coefficients of the master integrals can be simply extracted by 
solving a system of numerical equations, rather than computing phase-space integrals. Cut-by-cut, 
in a top-down algorithm, from quadruple- to single-cut, one can establish a system for extracting 
the coefficient of each master-integral identified by the product of the cut-propagators. The general 
structure of the integrand f3(| determines the polynomial shape of the equations forming such a sys- 



tem. In fact, by decomposing the loop variable in terms of a suitable basis of momenta (constructed 
from the external momenta and arbitrarily chosen reference momenta), the cut-conditions impose 
kinematical bounds on the values of the components of the loop momentum in this basis. According 
to the number of cuts, some component (when not all), are completely frozen, while others remain 
as free- variables 1 . The integrand, evaluated at a value of the loop momentum chosen among the 
solutions of a given cut, is a polynomial whose variables are the components of the loop momentum 
not frozen by the cut-conditions: the zero-th order coefficient of such a polynomial corresponds 
to the coefficient of the master integral one is interested in; the other terms of the polynomial, 
addressed to as spurious terms, do not contribute to the cut (since they vanish upon integration), 
but the are needed later on for extracting the coefficients of lower-point master integrals. 
To extract all the coefficients of this polynomial, one establishes a system, as said above, generated 
by evaluating numerically the integrand for values of loop momentum chosen within the solutions 
of the cut-conditions (parametrized by its free-components). The number of numerical evaluations 
must be the same as the number of the unknown coefficients to be determined. 

The goal of the current paper is to exploit the polynomial structure of the integrand and the 
freedom in choosing the solutions of the cuts, used as numerical points for the evaluation of the 
integrand, to improve the system-solving algorithm. By selecting the variables of each polynomial 
to be proportional to the roots of unity, the extraction of the polynomial's coefficients is carried 
through projections, using the orthogonality relation among wave-planes, rather than by system 
inversion. The basic principle underlying this procedure is the same as for the Discrete Fourier 



lr The variables not frozen by the cut-conditions correspond to the integration variables of the phase-space integral 
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Transform. The solutions accordingly obtained may help in getting a substantial gain in computing- 
time. The effective benefit of the new implementation is mostly experienced for the extraction of 
the coefficients of the 3- and 2-point functions, since in the 4- and 1-point cases the polynomial 
structure is very simple: a degree-1 polynomial, in the case of the 4-point; an effective degree-0 
polynomial, in the case of the 1-point functions. 

The paper is organized as follows. In Section 2, we recall the basic features of the OPP- 
reduction, by classifying the polynomial structures generated by the the multiple-cut integrand. 
In Section 3, we introduce the projections used for extracting the coefficients out of a polynomial, 
which are explicitly applied, in Sections 4 and 5, for the extraction of the coefficients of the 3- 
and 2-point functions, respectively. Finally, in Section 6, we apply the optimized reduction to 
the computation of the Next-to-Leading-Order QCD corrections to the scattering amplitude of the 
process ud — > W + W~W + . 



2. OPP-Reduction 



The starting point of the OPP reduction method |L8[ [19| is the general expression for the integrand of 
a generic m-point one-loop (sub-)amplitude that, using dimensional regularization, can be written 

as 

N(q) 



Di = (q + Pi) - m i , po^O. 



(2.1) 



D D 1 ---D m _ 1 

We use a bar to denote objects living in n = 4 + e dimensions; therefore we have q 2 = q 2 + q 2 , 
where q 2 is e-dimensional and (q ■ q) = 0. N(q) is the 4-dimensional part of the numerator of the 
amplitude. If needed, the e-dimensional part of the numerator should be treated separately, as 
explained in |35], |2l]]. N(q) depends on the 4-dimensional denominators A = (q + Pi) 2 — mf as 
follows 

m—1 m—1 

N{q) = ^2 [^(^0*1*2*3) + d(q; 20*1*2*3)] Y\. Di 



20<*1<*2<*3 
m—1 



j^io, 11,12,13 



m—1 



+ [cikhiz) + c(q; i iii 2 )] \ { A 



J0<»1<«2 



m—1 



m—1 



+ j2 [k*o*i)+%*o*i)] n a 

m—1 m—1 

+ ^2 K^o) + %; «o)] n A 
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m—1 



+ P(q)l[D l . (2.2) 

i 

Inserted back in Eq. ( |2.1| ), this expression simply states the multi-pole nature of any m-point one- 
loop amplitude, that, clearly, contains a pole for any propagator in the loop, thus one has terms 
ranging from 1 to m poles. 

The coefficients of the poles can be further split in two pieces. A piece that still depend on 
q (the terms d, c, b, a) , that vanishes upon integration, and a piece that do not depend on q (the 



terms d,c,b,a). Such a separation is always possible, as shown in Ul8| , and, with this choice, the 
latter set of coefficients is therefore immediately interpretable as the ensemble of the coefficients of 
all possible 4, 3, 2, 1-point one-loop functions contributing to the amplitude. Notice that the term 
with no poles, namely that one proportional to P{q) is polynomial and vanishes upon integration in 
dimensional regularization. Moreover, it can be shown that in the renormalizable gauge P{q) = 0, 
even before integration. 

2.1 Top-Down Polynomial Structures 

Since the scalar 1-, 2-, 3-, 4-point functions are known, the problem of computing the one-loop 
amplitude is simply reduced to the algebraical problem of fitting the coefficients d, c, b, a by eval- 
uating the function N(q) a sufficient number of times, at different values of q, and then inverting 
the system. 

This task can be achieved very efficiently by singling out particular choices of q such that, 
systematically, 4, 3, 2 or 1, among all possible denominators Di vanish. In ]T8 |, it was shown that 



by proceeding top-down from quadruple-cuts to single-cuts, it is possible to construct a particularly 
simple system of equations that can be solved analytically, whose solutions yield the complete 
reconstruction of the unknown coefficients. 

• Quadruple-cut. When q is solution of 

D = D 1 = D 2 = D 3 = 0. (2.3) 



N(q) = [d(0123) + d(q; 0123)] JJ A(<?) 

i^0,l,2,3 

= R{q) J] D t (q) (2.4) 

i^0,l,2,3 



where R(q) has a polynomial structure with 2 terms. 
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Triple-cut. At this stage all d and d coefficients are known. When q is solution of 

D = Di = D 2 = and D< ^ Vi ^ 0, 1, 2 (2.5) 

Eq. reads 

- ^K 012i 3) + %;oi2z 3 )] J] 

2<i 3 i^0,l,2,i 3 

= #(g) H A(g) = [c(012) + c(g;012)] J] , (2.6) 

1^0,1,2 #0,1,2 

where R'(q) has a polynomial structure with 7 terms. 

Double-cut. At this stage all d, d, c and c coefficients are known. When q is solution of 

A = Di = and A ^ Vz ^ 0, 1 (2.7) 

Eq. Q reads 

iV(g) - ^ [d(01 z 2 z 3 ) + Of i 2 t 3 )] J] ^ 

1<«2<«3 1^0,1,12,13 

- £[c(01i 2 )+c(g;01i 2 )] J] A 

1<«2 #0,1,12 

= i?"(g) J] A(g) = [6(01) + 6(g; Of)] J] A(<?) , (2.8) 

#0,1 #0,1 

where R"{q) has a polynomial structure with 9 terms. 

Single-cut. In massless theories all f -point functions, namely all tadpoles, vanish, also 
implying that, in such cases, one does not need to know all the b coefficients. However, in 
general, also the coefficients of the tadpoles are required. Therefore we discuss how to extract 
them. 

At this stage we assume to know all the d, d, c, c, b and b coefficients and, when q is solution 
of 

A = and A ^ Vz ^ , (2.9) 
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N(q) - Yl W^h) + d(q; Ozi^s)] \{ A 

0<il<i2<«3 1^0,11,12,13 



[c(0iii 2 ) + c(g; Oi!Z 2 )] ] | A 




5>(0ii) + 6(?;0ii)] |~] A 




iT (g) J] = [a(0) + %; 0)] J] A(?) • 



(2.10) 




where R"'(q) has a polynomial structure with 5 terms. We notice that the spurious coefficients 
aj(0) (i = 1,...,4) are never needed, because they would be necessary only to extract what 



extraction of the tadpole-coefficient a(0). 

We remark, that the polynomials we are going to deal with, i.e. R(q) , R' (?) , R" (q) , and R"'(q), 
share a common structure: a spurious structure, depending on the loop variable, q, which does not 
contribute to the cut-integral; and a single g-independent term, which corresponds to the actual 
coefficient of the master integral identified by the cuts 2 . 

3. Polynomial Structures and Discrete Fourier Transform 

After the general structure of Eq. ( |2.2|) is established, as we illustrated in the previous section, 
the calculation of the scattering amplitude reduces to the problem of extracting the coefficients of 
multivariable polynomials, generated at every step of the multiple-cut analysis. 

Let us show how it is possible to extract efficiently the coefficients of a polynomial of degree n 
in the variable x, say P n (x), defined as, 



by means of projections, according to the same principle underlying the Discrete Fourier Transform 
which works as follows. 

2 The only exception to this pattern will be the structure of the double-cut coefficients, that also includes a 
g-independent non vanishing term, in order to avoid numerical instabilities [19]. 



we called P(q), that, as already observed, is irrelevant. Therefore one focuses directly on the 




(3.1) 
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3.1 Discrete Fourier Transform 

Consider the function F, known only numerically in N points, F n (n = 0, N — 1). Each of this 
values, admits a DFT, defined as, 

N-l 

F n = J2h e~ 2m * n ■ (3.2) 

k=0 

The coefficients can be found by using the orthogonality relation, 

JV-l 

£ e 2m ^ n e~ 27Tl ^ n = N 8 ky . (3.3) 

n=0 

with the result 

N-l 

fk = ^Y< F " e2mJk 11 ■ ( 3 - 4 ) 

ra=0 

3.2 Projections 

Let us see how the above procedure can be implemented for extracting the coefficients q's of the 
polynomial in Eq.( |3.1| ) by projections. 

1. Generate the set of discrete values g k (k = 0, n), 

n 

9k = P n (x k ) = X> ^ e" 27ri (^) £ , (3.5) 

£=0 

by evaluating P n {x) at the points 

Xk = P e~ 2m ( n +v . (3.6) 

2. Using the orthogonality relations for the wave planes, one can obtain the coefficient q simply 
by, 



P 



-£ n 

-J^dk e 2m ^ £ (3.7) 



n + 

fc=0 



This procedure is very general and can be applied as long as one needs to know the coefficients 
of any polynomial. In fact, it has been recently suggested in ||29|| , for computing the coefficients of 
polynomials in the (— 2e)-dimensional mass parameter. 



S 



The projections could be extended also to the case of multi-variables polynomials, along the same 
line of the multi-dimensional DFT. Since we aim to minimize the computational load, we keep the 
number of numerical evaluations of each polynomial to be the same as the number of its coefficients. 
To this aim, we will see that the coefficients of a multi- variable polynomial can be equivalently found 
by breaking it in several one- variable polynomials, obtained from the former by freezing the values 
of the other variables, yielding still the use of the (one-dimensional) projections described above. 

The next two sections will be devoted to the application of the projection procedure for the 
extraction of the coefficients of the 3- and 2-point functions, respectively out of R'(q), and R"(q). We 
won't discuss hereby the reconstruction of the coefficients of the 4- and 1-point functions, because the 
polynomial structures of R(q) and R"'(q) is very simple, and the result of the projection procedure 
would be the same as the one given in ]18| . 



4. The coefficient of the 3-point functions 

In this section, we show how to apply the properties of orthogonal functions for extracting the 
coefficient of the 3-point functions. 



Let's begin from Eq.(2.21) of |8 



R'(q) = c(012) + {cy(012)[( g + po) ■ £ 3 } j + c 2j (012)[(q + p ) ■ i^} (4.1) 

3=1 



where R'(q) appeared in Eq.( [2.6| ). By substituting the parametrization of q, solution of the triple-cut 
conditions given in fl8|| , and recalled in App. 



q = —Po + %ih + + x 3 £ 3 + x 4 £ 4 (4.2) 

one obtains 

3 

R'{q) = c(012) + (M 012 ) & • Q J 4 + <%(012) (£ 3 • 4} , (4-3) 

3=1 

This expression can be read as a polynomial in x 3 and X4, whose canonical form reads, 

P(x 3 , x 4 ) = c ,o + ci )0 x 3 + c 2 ,o x\ + c 3i0 x\ + c ,i x 3 + c ,2 x\ + c ,3 x\ (4.4) 

with the following relations among the coefficients 

c ,o = c(012) , 

c jfi = c 2i (012) (4- Q j 

coj = c y (012) (£ 3 -£ 4 y (4.5) 
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Therefore computing the OPP-coefficients is equivalent to the computation of the 7 coefficients of 
Eq.(fyj). 



4.1 Projections 

The extraction of the coefficients Cjj is performed, in the framework of the original OPP-method and 
in CutTools, by choosing a redundand set of solutions in order to avoid some fake singularities oc- 
curring in kinematical points in which C = 0. In practice, this is obtained by doubling the number 
of calls to the numerator function N(q), by roughly doubling the computation time. 
The same problem can be more efficiently solved with the help of the proposed projection method. 
By using it, it is in fact very easy to find a set of solutions for which the point C = is never 
singular and for which N(q) is called just as many time as the number of needed coefficients. This 
is done at the price of creating a new fake singularity at C = 1, but, when this situation occurr, 
one can use the original solution. The described procedure explicitly shows the flexibility of the 
proposed projection method. 

To extract the 7 coefficients c^- by projections, we take 7 values of P, grouped in 2 sets (4 + 3), 



9i,k = P(x 3k ,x 4 k) , 

9l,k = P{x3k, ^4fc) , 



x 3k = C e" 2 ^ , x 4k = e 2 ^ , (k 



0,1,2,3) 
0,1,2) 



(4.6) 



where C is given in App. [A]. Then, we construct the auxiliary functions 




(4.7) 



fc=0 



with: (m,n) = (0,0), (0,1), (0,2), (1,0). 




(4.8) 



fc=0 



with: (m,n) = (0,0), (0,1), (1,0). 

In terms of these auxiliary functions, the coefficients read, 



Co,o = 0, 0) 



(4.9) 



1 



- /i(l, 1, 0)C n - 0, 0)C 8 + n(2, 0, 0)C 8 - /i(l, 0, 1)C 5 + /i(2, 0, 1)C 4 



Cl,0 = 




(4.10) 
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c 2 ,o = -^ty ( " Ml, 0, 2)C 10 + /i(2, 1, 0)C 8 - //(l, 1, 0)C 7 - /i(l, 0, 0)C7 4 

+/i(2, 0, 0)C 4 - Ml, 0, 1)C + /i(2, 0, 1) ) (4.11) 
c 3 ,o = —^—y ( " Ml. 0, 1)C 9 + ^(2, 0, 1)C7 8 - /i(l, 0, 2)C 6 + n(2, 1, 0)C 4 

-/i(l, 1, 0)C 3 - /i(l, 0, 0) + M2, 0, 0) ) (4.12) 
c 0) i = -7^—7 ( - /i(2, 0, l)C n + //(I, 0, 2)C 9 - fi(2, 1, 0)C 7 + /i(l, 1, 0)C 6 + /i(l, 0, 0)C 3 



C 12 _ I 



/i(2,0,l)C 2 + Ml,0,2) (4.14) 



-M2,0,0)C d + Ml,0,l)J (4-13) 
c ,2 = ( " M 2 > X > O)^ + Ml, 1, 0)C 9 + Ml, 0, 0)C 6 - /i(2, 0, 0)C 6 + //(l, 0, 1)C 3 



c ,3 = -^yU, 0, 0)C 9 - /i(2, 0, 0)C 9 + Ml, 0, 1)C 6 - /i(2, 0, 1)C 5 + /i(l, 0, 2)C 3 

-M2,1,0)C + //(1,1,0)J (4.15) 

Finally, to obtain the OPP-coefficients, simply use Eqs.( f4.5|) . 



5. The coefficient of the 2-point functions 

In this section, we show how to apply the properties of orthogonal functions for extracting the 
coefficient of the 2-point functions. 



Let's begin from Eq.(B.7) of [19 



R"(q) = b + b [(q + p )-v]+b 00 [(q + p )-v] 2 
+hi [(q + Po) ■ ii] + 621 [{q + Po) • 4] 
+M(g + Po) -ii} 2 + b 2 2[{q + Po) -is} 2 

+~boi [(q + Po) ■ £7} [{q + Po) ■ v] + ~b 02 [{q + Po) • £*) [(q + Po) ■ v) (5.1) 
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where R"(q) appeared in Eq. (|2.8|) . By substituting the parametrization of q, the solution of the 
double-cut conditions given in Eq.(B.6) of [19|], and recalled in App. ||, 



q = ~Po + yh + y v v + y 7 t 7 + x$£ s (5.2) 

one obtains 

N(q) = b + b [yh ■ v] + b m [yh ■ v} 2 

+fcii[y 8 4^7] + &2iM 7 -4] 

+~boi [yds ■ £7] [yh ■ v] + b 02 [y 7 t 7 ■ is] [yh ■ v] (5.3) 

This expression can be read as a polynomial in y, y 7 and ys, whose canonical form reads, 

P(y, y7, ys) = aooo + aioo y + a 2 oo y 2 

+«oio 1/7 + CJ020 y 7 + a oi ys + a 02 yl 

+«iio y 2/7 + ai M y y 8 , (5.4) 
with the following relations among the coefficients 

a ooo = b 

aioo = b (h • v) 

CI200 = b m (h ■ vf 

CtOlO = &2l(^7 • 4) 
O020 = &22^7 -4) 2 
«001 = &ll(4 -4) 

CJ002 = ^12(^7 -4) 2 

Olio = M^7 -4) {h •«) 

aioi = &oi(4 -4)(&i •«) (5.5) 

5.1 Projections 

Here we explicitly illustrate a solution that avoids the problem of doubling the number of calls to 
the numerator function N(q) (used in the original implementation of the OPP-method) due to the 
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appearence of fake singularities when F = 0. 

To extract the 9 coefficients by projections, we take 9 values of P, grouped in 5 sets (3 + 2 + 
2 + 1 + 1), 

g mh = P(0, y 7h , y 8h ) , y = 0, y 7h = F e" 2 ™t , y 8h = e 2 ™t , (h = 0, 1, 2) 
g 0h0 = P(0, y 7h , y Sh ) , y = , y 7h = e~ 2 ^ , y 8/l = F e 2m t (ft = , 1) 
0_ioa = P(-l, 2/7h, ys/,) , 2/ = "I , 2/7ft = P-i e" 2 ^ , y 8h = e 2 -f (ft = 0, 1) 
( 7 _ 1 = P(-1,1,F_ 1 ), 2/ = -l, 2/ 7 = l, j/ 8 = F_j 
( 7l = P(l,F 1 ,l), y = l , y 7 = F x , 2/ 8 = l 



where the definition of F y is given in App. 
Then, we construct the auxiliary functions 



i 2 

0, m, n) = - 9ooh e 27ri >- n > , (m, n) = (0, 0), (0, 1), (1, 0) (5.7) 



h=0 



2 

h=0 



2 

/i=0 



(5.6) 



/x(2, 0, m, n) = - £ g oh0 e 2 " 4 *^ , (m, n) = (0, 0), (1, 0) (5.8) 



/i(l, -1, m,n) = -J2 9-iok e 2 ^t(— «) , ( m , n ) = ( , 0), (0, 1) (5.9) 



M2, -1,0,0) = ^! (5.10) 



Ml, 1,0,0)= 9l (5.11) 
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The coefficients can be expressed as linear combinations of these auxiliary functions. Since the 
expressions for generic values of F Q , F_i, F\ are rather long, we present the one obtained when 
F = F_i = 0, which happens when the propagators are massless, as in the case of the application 
later discussed. Hence, the coefficient read, 



0000 =^(1,0,0,0) (5.12) 

1 

2 



aioo = ~ ( Ml, 0, 0, 0)F 2 - M2, 0, 0, 0)F 2 - Ml, -1, 0, 0)F 1 + /i(l, 0, 0, 0)Fi 



+M1, 0, 1, 0)F 1 + fi(2, -1, 0, 0)F 1 - /i(2, 0, 0, 0)^ - 2//(2, 0, 1, 0)F 1 

-Ml, -1, 0, 0) + Ml, -1, 0, 1) - 2M1, 0, 0, 1) + Ml, 1, 0, 0)) , (5.13) 



a 200 = I ( Ml, 0, 0, 0)F 2 - M2, 0, 0, 0)F 2 - //(l, -1, 0, 0)F 1 + Ml, 0, 0, 0)F, 



2 



+/i(l, 0, 1, 0)F 1 + /i(2, -1, 0, 0)Fi - //(2, 0, 0, 0)^ - 2/i(2, 0, 1, 0)^ 
+/i(l, -1, 0, 0) + Ml, -1, 0, 1) - 2/i(l, 0, 0, 0) 

-2M1, 0, 0, 1) - 2M1, 0, 1, 0) + Ml, 1, 0, 0) ) , (5.14) 



a io =M2, 0,1,0), (5.15) 

a 020 = /i(2, 0, 0, 0) - Ml, 0, 0, 0), (5.16) 

a 00 i = Ml, 0,0,1), (5.17) 

a 002 = Ml, 0,1,0), (5.18) 

a 110 = Ml, -1, 0, 0) - Ml, 0, 0, 0) - Ml, 0, 1, 0) - M2, -1, 0, 0) + M2, 0, 0, 0) + M2, 0, 1, 0), (5.19) 

o 101 - Ml, 0,0,1) -Ml, -1,0,1) . (5.20) 

Finally, to obtain the OPP-coefficients, simply use Eqs.( |5.5|) 



6. An example: QCD virtual corrections to ud — > W + W W + at NLO 

As an example of application of the optimized algorithm, we present the results for the pentagon 
and box diagrams contributing to the scattering amplitudes of ud — > W + W~W + at NLO in QCD. 
The complete cross section for the production of W + W~W + at LHC, as well as W + W~Z, W + ZZ 



and ZZZ is presented in [34 
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The main purpose of this application is to test the improvements on the OPP-reduction algo- 
rithm, both in terms of stability and efficiency, after the changes proposed in this paper have been 
included. 

The complete calculation of NLO QCD virtual correction to ud — > W + W~W + , neglecting the 
contributions that depend on the Higgs boson, involves the reduction of 53 diagrams. The topologies 
of pentagon and box diagrams contributing to this process are depicted in Fig. |1|. Overall, we have 
2 pentagons and 12 boxes. 




Figure 1: Box and pentagon diagrams contributing to virtual QCD corrections to ud — > W + W~W + . 
Wavy lines can be either Z or photons, dashed lines represent W + and W~. Diagrams involving an 
exchange in the two final W + should also be considered. 

Most of the computing time in the calculation is spent in the evaluation of the coefficients 
for 3-point and 2-point scalar function arising from the reduction of the diagrams illustrated. For 
example, the reduction of the each pentagon implies the evaluation of ten sets of c-type coefficients 
and ten sets of b-type coefficients, involving each seven and nine coefficients respectively. It is 
therefore very important to have efficient routines to achieve this task. As a comparison, the 
evaluation of the d-type coefficients is only performed five times for each pentagon and involves 
only two coefficients. 

For the purpose of this test, we fix the external momenta to the specific phase-space configu- 
ration of Eq. (|6.1| ) and the polarization vectors for the vector bosons. To perform the calculation, 
we use the OPP-reduction method, with and without the optimization technique illustrated before. 
The coefficients determined in this manner should be multiplied by the corresponding scalar inte- 
grals. Since, in the process that we are studying, no g-dependent massive propagator appear, we 
will only need massless scalar integrals. They are computed using the package OneLOop written by 
A. van Hameren J36|. Finally, we sum the contributions coming from the various pentagons and 



boxes: the results obtained are presented in Table ||. 

Pl = {500,0,500,0} 
p 2 = {500,0,-500,0} 

p 3 = {276.212,97.7237,-56.2856,238.9792} 
p A = {486.8926,213.4030,37.7214,-428.5282} 
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P5 = Pi + P2 ~ P3 ~ Pa 



(6.1) 



The results obtained with the two different implementations of the algorithm are, of course, in 
perfect agreement. However, with the optimized version we can roughly improve the efficiency of 
factor 2. As a first test, we checked the improvements on the timing of the two new routines alone, 
separating them from the rest of the reduction. We experience a reduction in the computing time 
of about 60% for the system of the c-coefficients, and about 50% for the system of the b-coefficient. 

In the overall evaluation of the amplitudes, the optimized routines are combined with other 
parts of the program that remain unchanged. This involves, for example, the initialization, the 
evaluation of scalar integrals, the evaluation of d- and a-coefficients and rational parts, and the 
summing over all contributions. We still retain, however, an overall improvement in the computing 
time of about 40%. 

Concerning the stability issues, we tested the new routines for a wide set of phase-space points. 
We do not observe significant improvements respect to the previous implementation. 



Polarization 


\A\ 





28.2435 


0- 


5.13851 


+ 


10.7870 


0-0 


1.923741 


- + 


0.718415 


+ 


7.43599 


+ - 


1.95506 


- + - 


0.276058 


- + + 


0.402302 


+ + + 


0.875457 



Table 1: Absolute value of the amplitudes for different configurations of the polarization of the vector 
bosons (accounting for 5-point and 4-point Feynman diagrams only). The results are expressed in unit of 
e 3 . 



7. Conclusions 

We presented an optimization of the reduction algorithm of onedoop amplitudes in terms of master 
integrals. That is based on the exploitation of the polynomial structure of the integrand when 
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evaluated at values of the loop-momentum fulfilling multiple cut-conditions, as emerged in the OPP- 
method. Accordingly, the integrand, evaluated at a value of the loop momentum chosen among the 
solutions of a given cut, is a polynomial whose variables are the components of the loop momentum 
not frozen by the cut-conditions: the zero-th order coefficient of such a polynomial corresponds to 
the coefficient of the master integral one is interested in; the other terms of the polynomial, though 
not contributing to the cut, are needed for the later determination of the coefficients of lower-point 
master integrals. 

To extract all the polynomial coefficients, one establishes a system generated by evaluating 
numerically the integrand for values of loop variable chosen within the solutions of the cut-conditions 
(parametrized by its free-components). The freedom in choosing the solutions of the cuts has 
been hereby exploited to improve the system-solving algorithm. By selecting the variables of each 
polynomial to be proportional to the primitive roots of unity, the extraction of the polynomial's 
coefficients is carried through projections, using the same orthogonality relation underlying the 
Discrete Fourier Transform. The number of numerical evaluations is kept as low as the number 
of the unknown coefficients to be determined, by using one-dimensional projections also in case of 
polynomials in more than one- variable. 

The novel implementation was applied to the reduction of the 4- and 5-point one-loop Feynman 
diagrams contributing the NLO QCD corrections to ud — > W + W~W + , where we experienced a 
reduction of the computational load. 

The flexibility of the projection-procedure hereby presented extends its range of applicability 
to tackle the determination of the coefficients of polynomial structures wherever should this issue 
occur. Moreover, we finally remark that the parametrization of the free (integration) variables 
as complex unitary phases yields as well a very effective performance of Cauchy's residue theorem 
within the contexts of factorization- and unitarity-based methods, where the on-shellness properties 
are naturally captured by polar structures in complex phases. 
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A. The basis for the 3-point functions 

The loop momentum solution of the triple-cut, given in Eq.( [4.2| ), is expressed in terms of auxiliary 
vectors defined as follows. 

l\ and £ 2 are massless 4-vector satisfying the relations 

fc 1 = 4 + ai 4 ; k 2 = £ 2 + a 2 h, (A.l) 

with 

ki=pi-p . (A. 2) 

Furthermore, in spinorial notation, 

%=(h\-f\h\, ^=<WN with (£ 3 ■ Q = -m ■ £2) ■ (A.3) 



The solution to Eq. flA.lf) reads 

£ x = (3(ki — a x k 2 ) , £ 2 = (3(k 2 -a 2 ki), 

k 2 

(3 = 1/(1 - aia 2 ) , on = — , 

7 

7 = 2{h ■ £ 2 ) = {k x ■ k 2 ) ± VA , A = (h ■ k 2 f - k\k\ . (A.4) 
we decompose g M + Pq in the basis formed by £ x ,£ 2 ,£3, and £4, 

Q = —Po + x\h + x 2 £ 2 + x 3 £ 3 + x 4 £ 4 , (A. 5) 

which is solution of the triple-cut, 

D — D x — D 2 — . (A.6) 

Due to the above constrains, the coefficients of the loop decomposition must fulfill the following 
relations, 

£1 = —[d 2 — a 2 d x — d (l — a 2 )) , (A. 7) 

7 

x 2 — —\d\ — a\d 2 — c? (l — 011)) , (A. 8) 

7 

x 3 x 4 = C , (A. 9) 

where 

C = \(xi*2-^j , (A.10) 
di = m 2 -k 2 . (A. 11) 
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B. The basis for the 2-point functions 

First, we introduce a massless arbitrary 4-vector v, such that (v ■ ki) ^ 0, that we use to rewrite k\ 
in terms of two massless 4- vectors (we also take £ 2 = 0) 

k 1= £ + av, (B.l) 

giving 

k 2 

7 = 2(h -v) = 2(£-v) and a = -± . (B.2) 

7 

Then, we introduce two additional independent massless 4-vectors £j^ defined as 

% = (l\>f\v] , e 8 = (v\Y\i], (B.3) 

for which one finds 

(£ 7 -£s) = -2j, (B.4) 
and we decompose g M + p$ in the basis formed by ki, v, £ 7 and £ 8 

<? = -p£ + vK + + vA + > (b.5) 

that fulfill the double-cut requirement 

D = D 1 = 0. (B.6) 
For a q written as in Eq. ( |B.5| ) this implies the system 

yivs = F y 

di-do- 2yk\ 

Vv = , (B.7) 

7 

where 

F y = -i- (m 2 -y(d 1 - d ) + y 2 kj) . (B.8) 
We remark that when — 0, F y vanishes for y — — 1, 0. 
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